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SUMMARY 


This  report  discusses  the  initial  modeling  effort  to  develop  a 
tool  which  can  assist  in  the  development  of  depot  maintenance  work 
requirements  by  determining  the  inherent  reliability  of 
equipment.  A  Monte  Carlo  simulation  was  developed.  The  input  for 
the  simulation  includes  a  block  diagram,  a  time  to  failure 
distribution  or  failure  and  suspension  data  for  each  block  in  the 
diagram,  and  the  age  of  each  block.  The  primary  output  of  the 
simulation  is  a  numerical  time  to  failure  distribution  for  the 
system.  This  distribution  is  used  to  generate  a  graph  of  the  time 
to  failure  distribution,  a  histogram  of  the  ages  at  which  the 
system  fails,  a  graph  of  hazard  rates,  a  total  time  on  test  plot 
and  an  graph  of  the  expected  cost  of  various  replacement  policies 
if  cost  data  is  available. 

The  algorithms  and  data  structure  used  in  the  model  are  discussed 
as  well  as  the  potential  uses  of  the  model  and  its  performance. 
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1. 


Introduction 


a.  Background 

This  effort  was  initiated  as  the  result  of  a  proposal  made  by 
Gerald  Moeller,  while  a  LOGAMP  trainee  on  assignment  to  the 
Maintenance  Directorate  at  the  Headquarters  U.S.  Armament, 
Munitions  and  Chemical  Command.  The  original  proposal  was  to 
construct  a  model,  which  would  support  overhaul  decisions  at 
depot  shops.  After  a  series  of  discussions  between  the 
Maintenance  Directorate  and  the  Systems  Analysis  Office,  it  was 
decided  that  the  model  should  be  pursued  as  an  aid  to  depot 
maintenance  work  requirement  development,  rather  than  a  tool  for 
the  depot  shop  foreman. 

b.  Objective 

The  initial  objective  was  to  develop  a  proof  of  principal  model 
capable  of  generating  hazard  rates  and  assessing  a  system's 
reliability  once  repairs  were  completed.  This  model  would 
provide  a  tool  to  assist  in  the  identification  of  components, 
which  should  be  rebuilt,  replaced  or  left  alone,  when  the  system 
is  overhauled.  The  model  had  to  be  structured  in  a  flexible 
manner  so  it  could  readily  host  the  modeling  of  any  weapon 
system  in  a  manner  which  enables  combining  the  reliabilities  of 
the  various  components  or  sub-assembilies  in  the  peculiar  way 
which  reflects  that  end-item's  true  reliability.  The  model  also 
had  to  be  constructed  in  a  manner  which  would  allow  personnel 
having  limited  computer  experience  to  use  it  productively.  The 
objective  of  this  report  is  to  summarize  the  modeling  efforts 
to  date. 
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2.  Methodology 


a.  Selection  of  Computer  Language  and  Method 

The  model  was  coded  in  Turbo  Pascal  because  it  supports 
recursion,  dynamic  memory  allocation  and  the  graphic  routines 
required  by  the  algorithms  used  in  the  model. 

In  order  to  achieve  its  objectives,  the  model  would  have  to 
support  the  parsing  and  analysis  of  block  diagrams  to  estimate 
system  reliability. 

Two  modeling  methods  were  considered,  Monte  Carlo  simulation  and 
conventional  probability  analysis.  Monte  Carlo  simulation  was 
selected  because  of  its  flexibility  and  the  difficulty  in 
automating  the  process  of  converting  a  block  diagram  into 
equations  for  models  more  complex  than  a  series/parallel  system. 
The  simulation  Cl ike  all  models)  is  a  blend  of  algorithms  and 
data  structures,  which  are  explained  in  the  main  body  of  this 

report.  In  addition  to  the  simulation,  the  program  also 

contains  a  unit  to  analyze  raw  data  and  a  unit  to  produce 

graphic  output.  The  algorithms  used  in  these  units  and  the 

algorithms  used  to  generate  the  random  variables  used  in  the 

simulation  are  given  in  the  appendixes. 

b.  Simulation  Algorithm 

The  Monte  Carlo  method  involves  random  sampling  from  the 
conditional  time  to  failure  distributions  of  the  components 
making  up  the  block  diagram  and  using  these  times  to  determine 
the  time  at  which  the  system  fails.  The  procedure  is  repeated 
many  times  until  enough  failure  times  are  obtained  to  estimate 
the  system  reliability.  The  resulting  data  can  be  used  to 
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FIGURE  1.  FLOWCHART  FOR  SIMULATION 


estimate  parameters  for  candidate  time  to  failure  distributions 
for  the  system.  A  flowchart  of  the  logic  used  in  the  simulation 


is  provided  in  Figure  i.  In  addition,  the  procedure  is 
illustrated  by  the  following  example.  The  mission  reliability 
diagram  is  given  in  Figure  2.  The  mission  requirement  would  be 
that  equipment  A  and  either  equipment  Ci  or  C2  work,  or  that 
equipment  B1  and  Cl  work,  or  that  B2  and  C2  work  for  success. 
The  time  to  failure  distribution  and  age  of  each  component  is 

FIGURE  2.  MISSION  BLOCK  DIAGRAM 


NAME: 

B 1 

DIST: 

EC  60000) 

AGE: 

26000 

NAME: 

Cl 

DIST: 

VC  O.  6.  30000.  O) 

AGE: 

26000 

NAME: 

DIST:  UN  C  1  0 .  6,  1,03 

AGE:  26000 


NAME: 

B2 

DIST: 

E  C  60000 5 

AGE: 

26000 

NAJA: 

C2 

DIST: 

WC  1 . 6.  64000,  05 

AGE: 

26000 

provided  in  each  block.  The  age  of  components  can  be  measured 
in  any  appropriate  unit,  such  as  hours  or  miles,  but  the  units 
must  be  the  same  for  each  block.  Source  and  sink  blocks  have 
been  added  to  the  diagram  because  they  are  required  by  the 
model’s  data  structure.  The  source  block  starts  the  simulation 
and  the  sink  block  ends  the  simulation. 

The  model  first  generates  random  numbers  between  0.0  and  1.0  and 
uses  them  ,  with  the  algorithms  in  Appendix  B, .  to  calculate  a 
time  to  failure  for  each  component.  These  times  are  given  in 
Table  1. 

TABLE  1.  TIME  TO  FIRST  FAILURE  FOR  COMPONENTS. 


BLOCK 

TIME 

A 

3162. 36 

B1 

46487.67 

B2 

10291.19 

Cl 

103741.06 

C2 

97571.61 

4 


The  component  failure  times  are  sorted  and  the  system  time  set 
equal  to  the  minimum  component  failure  time.  Using  the 
FIGURE  3.  FLOWCHART  FOR  PROCEDURE  EVALUATE  GIVEN  BLOCK 


procedure  Evaluate  Csee  Figure  3.),  the  model  determines  if  a 
success  path  from  source  to  sink  can  be  found  among  the  failed 
and  nonf ailed  equipments.  The  model  first  assigns  a  value  of 
true  to  the  boolean  variable  SYSTEMFAILS.  It  then  evaluates  the 
blocks  by  comparing  their  time  to  fail  with  the  system  time.  If 
the  time  to  fail  is  greater  then  the  system  time,  the  procedure 
also  evaluates  the  blocks  for  which  the  current  block  is  the 
origin  node.  If  the  block  being  evaluated  is  the  sink  block  and 


3 


it's  time  to  failure 

is 

greater 

than 

the 

system 

time. 

the 

variable  SYSTEMFAILS 

is 

set 

to 

false. 

The 

system 

time 

is 

updated  to  the  time 

of 

the 

next 

component 

failure 

and 

the 

process  continued  until  a  success  path  does  not  exist.  Table  2 
shows  the  various  system  times  and  the  condition  of  the  system 
at  each  time. 


TABLE  2.  SUCCESS/FAILURE  ARRAY  FOR  BLOCK  DIAGRAM 


SYSTEM  TIME 

A 

Bl_ 

B2 

Cl 

C2 

SYSTEM 

3162.36 

F 

S 

S 

s 

S 

S 

10291.19 

F 

s 

F 

s 

s 

S 

46487 . 67  * 

F 

F 

F 

s 

s 

F 

97571.16 

F 

F 

F 

s 

F 

F 

105741.06 

F 

F 

F 

F 

F 

F 

*  The  system 

fails 

at 

this 

point. 

No 

futher 

evaluation 

would 

take  place  during 

this 

iteration. 

S  = 

SUCCESS 

F  =  FAILS 

When  the  system  fails,  the  model  records  the  age  of  the  system 
at  failure,  updates  counters  and  repeats  the  process  until  the 
number  of  iterations  performed  equals  the  number  required.  Once 
the  required  iterations  are  completed,  the  ages  are  sorted  and 
the  associated  probabilities  estimated,  using  median  ranks. 
The  two  keys  to  performing  this  simulation  are  the  recursive 
procedure  Evaluate  and  the  data  structure,  which  allows  the 
model  to  traverse  the  block  diagram. 

c.  Data  Structure 

The  pri»  . y  building  blocks  of  the  simulation's  data  structure 
are  record'  jsociated  with  each  component.  These  records  are  a 
combinat Ion  of  simple-type  data  into  a  new  data  type,  which  has 
two  advantages.  All  data  elements  for  a  single  record  are 
logically  connected  to  each  other.  Also,  some  operations,  such 
as  assignment,  can  be  performed  on  the  entire  record, 
eliminating  the  need  to  refer  to  each  element  of  the  record. 
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These  records  are  made  up  of  two  kinds  of  data  elements,  those 
that  contain  information  about  a  system  component  and  those  that 
point  to  other  records.  The  pointers  connect  blocks  to  form  a 
linked  list  for  data  editing  and  a  tree,  which  is  used  to 

traverse  the  block  diagram.  Appendix  D  contains  the  data 
elements  which  make  up  each  record  and  a  brief 

discussion  of  how  they  are  used. 

3.  Data  Requirements 

The  data  requirements  for  this  or  any  model  meeting  the  above 

objectives  include  the  following. 

a.  A  block  diagram  for  the  equipment  to  be  overhauled. 

b.  A  time  to  failure  distribution  for  each  block  or  failure  and 

suspension  data  which  can  be  used  to  generate  the  failure 

distributions.  The  failure  distributions  currently  available 
are  exponential,  uniform,  weibull,  normal,  lognormal  and  a 
constant  probability  of  failure. 

c.  The  age  of  each  block. 

d.  The  cost  of  scheduled  replacement  and  of  field  failures  if 
an  expected  cost  distribution  is  required. 

Availability  of  this  data  is  the  major  obstacle  to  the 

implementation  of  this  type  of  model. 

4.  Model  Performance 

a.  Accuracy 

All  Monte  Carlo  estimates  have  an  associated  error  band.  The 
larger  the  number  of  iterations  the  more  precise  the  estimate 
is.  A  100C1.0  -  c*>  percent  confidence  interval  <CI>  on  a 

given  reliability  estimate  can  be  constructed,  using  the 
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following  formula,  which  comes  from  the  normal  approximation  of 
the  binomial  probability  density  function. 


/*  ”  * 
p  Ci  -  p> 

N 

where  '2.<xyz  is  the  100(1  -  a/2)  percentile  point  of  the  standard 
normal  distribution,  p  is  a  point  estimate  of  the  reliability 
and  N  is  the  number  of  iterations. 

The  model’s  output  for  the  bridge  diagram  in  Figure  1.  was 
compared  with  the  actual  time  to  failure  distribution  and  the 
maximum  absolute  error  (MAC)  recorded.  The  results  are  given 
in  Table  3  with  the  largest  error  expected  for  90  and  95  percent 
confidence  intervals. 


Table  3.  MODEL  OUTPUT  VS  ACTUAL  RELIABILITIES 


— Baa iff On 

WIDTH - 

#  Iterations 

MAE 

90k  Cl 

95k  Cl 

1000 

0.02476 

0. 03010 

0.03542 

2000 

0.01505 

0.02191 

0.02504 

3000 

0.01472 

0.01789 

0.02045 

4000 

0.01356 

0.01550 

0.01771 

5000 

0.01625 

0.01386 

0.01584 

6000 

0.01222 

0.01265 

0.01446 

7000 

0.01144 

0.01171 

0.01339 

A  confidence  interval  on  the  error  in  the  estimate  of  the 

A 

time  to  failure  can  be  obtained  using  E  ■  ±  ZQ^2  —  , 

N 


mean 

where 


a  is  an  estimate  of  the  standard  deviation. 

Since  this  is  a  proof  of  principal  model,  the  number  of 
iterations  has  been  preset  to  1000.  If  this  model  becomes  a 
standing  model  it  should  be  modified  to  allow  the  analyst  to  set 
the  number  of  iterations  as  well  as  the  random  number  seeds  and 
other  parameters  making  up  the  simulation  environment.  Also 
variance  reduction  techniques,  such  as  stratified  sampling, 
should  be  employed. 


b.  Execution  Time 
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The  length  of  time  it  takes  to  run  the  model  depends  on  the 
number  of  blocks  in  the  diagram,  the  distributions  used  and  the 


complexity  if  the  block  diagram.  The  run  times  for  various 
configurations  is  given  in  Table  4. 


Table  4  EXECUTION  TIME 


Configuration 

Time 

1. 

One  Block  with  various  Distributions 

4-6 

sec 

2. 

3. 

The  Bridge  diagram  in  Figure  2. 

Five  Exponential  Distributions  in 

22 

sec 

4. 

parallel 

Five  Exponential  Distributions  in 

19 

sec 

5. 

series 

100  blocks  in  series  with  Uniform 

distributions.  The  last  blocks  had 

a  maximum  time  to  failure  less  than 

14 

sec 

6. 

the  mimimum  of  any  other  block. 

The  same  conf iguratioin  as  number  5, 

3  min  93 

sec 

but  with  1000  blocks.  42  min  14 

sec 

5.  Potential  Uses  of  the  Model 

a.  Support  of  Overhaul  Decision  Making 

The  current  model  meets  the  initial  objectives  of  generating 
hazard  rates  and  assessing  a  system's  reliability  once  repairs 
are  completed.  It  can  be  used  to  support  decision  making 
concerning  when  Cor  if)  overhaul  should  take  place,  by 
generating  time  to  failure  distributions  for  various  levels  of 
overhaul.  The  model  can  support  decisions  based  on  mission 
reliability,  change  in  the  MTTF  or  on  the  basis  of  cost. 
However,  the  current  configuration  of  the  model  requires  several 
runs,  with  the  analyst  manually  modifying  the  data.  A  possible 
enhancement  to  the  model  would  be  the  generation  of  a  list  of 
components  ordered  by  the  effect  of  their  overhaul  on  the  the 
system  reliability. 
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b.  General  Tool  for  Reliability  Analysis. 


Because  the  information  generated  by  the  model  is  fundamental  to 
any  analysis  of  the  reliability  of  a  system  or  component,  the 
model  can  be  used  to  support  reliability  analysis  of  block 
diagrams  or  failure  and  suspension  data.  The  model  can  generate 
the  following  information  for  a  block  diagram. 

(1)  A  numerical  x.ime  to  failure  distribution 

<2)  Numerical  conditional  time  to  failure  distribution 

(3)  Numerical  reliability  distribution 

(4)  Hazard  rates 

<5>  Mean  time  to  failure 

(6)  The  probability  each  component  outlives  the  system 
The  model  can  support  the  following  data  analysis. 

(1>  Estimation  of  parameters  for  time  to  failure 

distributions. 

<2>  Nonparametric  estimation  of  hazard  rates. 

(3)  Calculation  of  mortality  curves  using  ranking. 

c.  A  Base  for  Further  Modeling 

The  current  model  provides  a  foundation  for  models  intermixing 
maintenance  policies  and  the  inherent  reliability  of  hardware. 
Such  models  could  be  constructed  in  two  ways.  The  preferable 
method  would  be  to  use  the  data  generated  by  the  model  to 
construct  a  time  to  failure  function.  This  could  be  done  using 
the  current  raw  data  analysis  unit  of  the  model  or  by  fitting 
various  functic  is  to  portions  of  the  data  and  generating  a 
function  in  a  piecewise  manner.  Once  a  distribution  was 
constructed,  it  could  be  used  (with  other  information)  to 
generate  such  things  as  the  effect  of  corrective  and  preventive 
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maintenance  on  the  reliability  of  a  system  and  estimates  of  the 
steady  state  availability  of  maintained  systems. 

The  second  method  would  be  to  continue  to  build  on  the  existing 
simulation.  Maintenance  actions  would  be  scheduled,  in  the 
simulation,  in  a  fashion  similar  to  the  way  component  failures 
are  scheduled  in  the  current  simulation.  The  effect  of 
maintenance  actions  would  be  modeled  by  generating  new  failure 
times  for  the  components  on  which  maintenance  was  performed  and 
recording  such  things  as  down  time  to  allow  the  model  to 
estimate  availability  as  well  as  the  effect  of  the  maintenance 
on  reliability.  Such  a  model  would  allow  analysis  of  systems 
too  complex  to  work  with  analytically. 

6.  Conclusions 

The  development  of  the  proof  of  principal  model,  demonstrates  the 
analytical  feasibility  of  computer  models  performing  reliability 
analysis,  which  can  support  decision  making  regarding  the  type  of 
maintenance  that  should  be  performed.  However,  such  analysis 
should  only  be  carried  out  in  conjunction  with  an  engineer 
familiar  with  the  equipment  being  studied.  An  analysis  without 
this  expertise,  could  easily  result  in  incorrect  conculusions  and 
poor  maintenance  planning. 

The  question  of  the  availability  of  accurate  data  to  run  such 
models  seems  to  be  the  major  objection  to  their  use.  The  lack  of 
good  data  can  be  a  problem,  especially  early  in  the  life  of  a 
weapon  system.  However,  plans  are  being  made  and  must  be  made 
with  the  data  available.  The  important  question  is  whether  the 
model  will  or  will  not  make  better  use  of  the  available  data. 

Ve  believe  it  will.  The  primary  product  of  the  model,  a  numerical 
time  to  failure  function,  is  fundamental  to  the  analysis  of  the 
reliability  and  maintainability  of  any  system.  In  all  but  the 
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most  simple  cases  reliability  analysis  entails  rather  complicated 
mathematical  formulations  and  tedious  computations.  This  type  of 
model  provides  a  structured  method  for  carrying  out  such  analysis 
with  a  minimum  of  simlifying  assumptions,  while  providing  a  tool 
for  the  systematic  analysis  of  failure  data  and  updating 
assessments  of  the  reliability  of  components  as  well  as  the 
system. 
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APPENDIX  A  MORTALITY  CURVE  CONSTRUCTION 


A  mortality  curve  gives  values  for  the  proportion  of  a  population 
that  fail  CP  >  before  a  given  time.  The  model  calculates  the 
points  on  a  mortality  curve  using  ranking.  It  requires  a  data  set 
consisting  of  the  age  at  which  failure  or  suspension  takes  place 
CX  ),  the  number  of  items  exposed  in  the  interval  -  Xv-1  CE^>, 
and  the  number  of  items  which  failed  in  the  interval,  X^  - 
<F  )  for  v  -  1,  2,  3,  .  .  N.  The  following  algorithm  is  used. 

a.  Set  R  -  1.0  and  MAXE  -  MaxCE^  E2,  .  .  .  ,  EN> .  Order  the 
data  chronologically. 


b.  Set  R 

c .  Set  P 


RC1  - 


)  and  P 


1  -  R,  for  l  “1,  2, 


N. 


MAXE  P  -0.3 


for 


1,  2, 


MAXE  +0.4 

For  uncensored  data,  this  algorithm  reduces 
which  is  a  commonly  used  appromiximation  of 


to  P 
the 


N. 

t  -  0.3 

i.  *  N  4-  on  ’ 

median  rank  of 


the  ith  failure. 
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APPENDIX  B  DISTRIBUTIONS  SUPPORTED  BY  THE  MODEL 


1.  This  appendix  contains  information  concerning  the 
distributions  supported  by  the  model.  For  each  distribution,  the 
following  information  is  provided. 

a.  A  brief  discription  of  how  the  distribution  is  usually  used 
in  reliability  modeling. 

b.  The  conditional  time  to  failure  distribution. 

c.  The  algorithm  used  to  generated  random  variables  from  the 
conditional  time  to  failure  distribution. 

d.  The  algorithm  used  to  estimate  the  distribution’s  parameters 
given  a  set  of  failure  and  suspension  data. 

2.  Exponential  FCx)  *  1  -  e"x*- 

a.  The  exponential  distribution  is  the  most  widely  used  in 
reliability  .  It  can  be  used  to  model  random  component  failures 
and  has  a  constant  failure  fate  of  A. 

b.  The  probability  of  failure  at  time  Cx+s)  given  the  component 
has  survived  to  time  x  is  given  below. 

PCX  >  x+s  jX  >  x)  -  PCX  >  x  +  s)  ✓  PCX  >  x> 

-\Cx  +  s) 


PCX  >  x+s|X  >  x)  -  e“AS 
PCX  s  x+s |X  >  x)  «  i  -  e~Xs 

c.  Random  variables  for  the  conditional  exponential 
distribution  are  generated  using  the  inverse  transform  method. 

Cl)  Generate  U  -  UCO.l) 

C2)  Set  s  -  X  InCU) 

d  An  estimate  of  X  is  made  by  applying  ordinary  least  squares 
COLS)  to  the  linerized  distribution  and  then  using  the  derivative 
of  the  sum  of  squared  errors  CSSE)  to  search  in  the  neighbourhood 

A 

of  X  for  an  estimate  of  X  which  minimizes  the  SSE. 


Cl)  Set  Y 
for  i  "1,  2,  . 


-InCI . 0  -  PCX  £  x  ));  and  x  -  InCX  ) 

V  V  V 

N  Where  N  equals  the  number  of  observations. 


Generate  X  for  Y  ■  A  X  using  OLS.  Set  L  ■  0  and  U  ■  2  A. 

N 

C2)  Set  -  2  “2  Xie”WCi  {PCX  £  x.  )  -  1  +  e"XXij- 


C3)  If  -S?E  *  0.0  then  set  U 
d  X 


Otherwise,  set  L  ■  X. 
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Set  \  -  L  +  CU  -  L)  /  2.0. 

C4)  If  —  —  s  0.00001  stop.  Otherwise,  go  to  2. 

X 

CX  -  a) 

3.  Uniform  FCX)  -  - 

(b  ~  a> 

a.  This  distrubition  is  not  used  very  often  in  reliability.  It 
could  be  used  in  cases  where  the  only  thing  known  about  a 
component  is  its  maximum  life. 

b.  PCX  >  x+s |X  >  x)  -  PCX  >  x  +  s>  ✓  PCX  >  x> 


PCX  >  x+s  |X  >  x) 
PCX  £  X+s |X  >  x> 


c.  Random  variables  from  the  conditional  uniform 
distribution  are  generated  using  the  inverse  transform  method. 

Cl>  Generate  U  -  UC0,1) 

C2>  set  s  -  U  Cb  -  x) 

d.  Estimates  of  the  parameters  of  the  Uniform  distribution  are 
made  by  applying  OLS  to  failure  and  suspension  data. 

Cl)  Set  Y  -  PCX  s  x  ). 


C2>  Estimate  a  and  b  for  Y^  ■  ^  1  xi  -  ^  ^  using  OLS. 

4.  Veibull  FCX)  ■  1  -  e“C 

a.  The  weibull  distribution  is  the  most  general  one  used  to 
model  reliability.  It  can  be  used  to  model  initial  failures  or 
wearout  failures.  Its  failure  rate  depends  on  the  shape 
parameter,  ft. 

Cl)  If  ft  >  1  increasing  failure  rate. 

C2)  If  ft  <  1  decreasing  failure  rate. 

C3)  If  ft  ■  1  constant  failure  rate. 


b.  PCX  >  x+s  JX  >  x> 
PCX  >  x+s |X  >  x) 


PCX  >  x+s)  ✓  PCX  >  x) 
e~  C  C  x+s-y)  Srjl  P 
e-CC  x-y)^]^ 
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PCX  s  x+s jx  >  X>  -  i  -  +  icx-r>/niP 

c.  Random  variables  from  the  conditional  weibull  distribution 
are  generated  using  the  inverse  transform  method. 

Cl>  Generate  U  ^  UC0,1> 

C2>  Set  s  -  77£-lnCl-U)  +  CC +  y  -  x 

d.  Estimates  of  y,  77  and  ft  are  made  by  combining  a  golden 
section  search  and  OLS  to  obtain  estimates  which  minimize  the  SSE. 

Cl)  Set  B  ■  0.0.  Set  U  -  minCX^  X2,  .  .  .  XN>. 

If  minCXt,  X2,  .  .  .  XN)  >  0  then  set  T  ■  minCX1>  X2,  .  .  XN> . 

Otherwise  set  T  m  1.0.  Set  y%  ■  0.382  CU  -  B>  +  B. 

Set  X^  *  lnCXt  -  7^  and  Y^  ■  ln-^ln  ^  ■_  ^  ^  for 

v  *  1,  2,  . . . ,  N.  Generate  OLS  estimates  for  Yt *  ftX^+  ft  lnCrp 

and  SSE4. 

C2>  Set  y  -  0.618  CU  -  B>  +  B.  Set  X  ■  InCX  -  >  >  and 

'2  1  1  2 

Y^  *  ln-^ln  jr£"7T 3 j-J-  for  c  ■  1,  2,  .  ..,  N.  Generate  SSE2 

and  OLS  estimates  for  Yv  ■  (J  X  +  ft  InCrp. 

C3)  If  SSE.  >  SSE  1 
1  2 

then  set  B  *  y±;  7t  m  7Z>  7Z  *  0.618  CU-B>  +  B;  X^-CX^-?^)  and 

SSE  *  SSE  .  Generate  OLS  estimates  and  SSE, 

12  2 

Otherwise,  set  U  ■  y2  ;  7zm  7 xi  7t  ■  0.383  CU  -  B>  B; 

X  *  CX  -  >  >  and  SSE,"  SSE..  Generate  OLS  estimates  and  SSE.. 
11*121  1 

|C >  -  7> | 

C4>  If  I —  -  >  0.0001  go  to  3.  Otherwise,  stop. 


5.  Normal 


FCX)  - 


-Cx-/i>2/2o2 

e 


dx 


a.  No  closed  form  exists  for  this  distribution.  Its  failure 
rate  is  a  monotonically  increasing  function  of  X.  The  normal 
distribution  can  be  used  to  model  weareout  and  stress  failure 
where  the  random  variable  is  stress  rather  than  time. 

b.  There  is  no  closed  form  of  the  conditional  time  to  fail 
distribtion. 

c.  Random  variables  from  the  conditional  normal  distribution, 

PCX  £  x+s|X  >  x),  are  generated  using  a  standard  normal 
transformation  and  the  inverse  transform  method. 
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Cl)  Generate 

distribution. 

Y  -  PCX  £  x  ) 

from 

the 

standard 

normal 

C2) 

Generate 

U  —  UC0,1)  and  set 

U  «  u' 

Cl  - 

Y)  +  Y 

C3) 

Generate 

Z  associated  with  U 

from 

the 

standard 

normal 

distribution. 


C4)  set  s  •  +  o  Z  -  x 

d.  Estimates  for  the  parameters  of  the  Normal  distribution  are 
made  by  applying  OLS  to  the  standard  normal  distribution 
associated  with  failure  and  suspension  data. 

Cl)  Generate  Y  as  the  inverse  standard  normal  value 

V 

associated  with  PCX  1  X  )  for  i.  •  1,  2,  .  .  .  N. 

C2>  Use  OLS  to  estimate  parameters  for  ^  and 

use  these  parameters  to  estimate  and  a. 


6.  Lognormal 


x  a  tz  n 


exp 


a.  The  failure  rate  of  the 


-Cln  x  -  p) 

To 5 

lognormal 


2 

—  dx 

distribution 


initially 


Increases  over  time  and  then  decreases,  approaching  zero.  It  is 


useful  for  modeling  situations  in  which  early  failures  dominate. 

b.  .  There  is  no  closed  form  of  the  conditional  time  to  fail 
distribution. 

c.  Random  variables  from  the  conditional  lognormal 

distribution,  PCX  >  x+s|X  >  x),  are  generated  using  an 

acceptance/rejection  method. 

Cl >  Generate  Y  «/«  N  C/u,o|X  >  x) 

C2)  Set  s  -  eY+x  -  x 

d.  Estimate  for  the  paramters  of  the  Lognormal  distribution  are 
made  by  estimating  the  parameters  of  the  Normal  distribution 
associated  with  the  Lognormal. 

Cl>  Set  X  -  InCX  )  for  v  -  1,  2,  ...  N. 

VI 

C2)  Generate  Y  as  the  inverse  standard  normal 

l 

value  associated  with  PCX  £  Xt )  for  v  ■  1,  2,  .  .  .  N. 

C3>  Use  OLS  to  estimate  parameters  for  Yt  ■  Xt  -  and 

use  these  parameters  to  estimate  n  and  o. 

7.  In  addition  to  the  above  distribution  the  model  also  supports 
use  of  components  with  constant  Cnot  age  dependent)  probabilities 
of  failure.  Random  variables  are  generated  for  the  components 
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using 

a. 

b. 

set  s 


appendix  b  distributions  supported  by  the  model 

the  inverse  transform  method. 

Generate  U  UC0,1> 

If  U  >  PC component  fails)  then  set  s  -  103O°-  Otherwise, 

-  0.0. 
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APPENDIX  C  GRAPHICS  GENERATED  BY  THE  MODEL 


FIGURE  4.  TIME  TO  FAILURE  DISTRIBUTION 


0.03  1.05  2.07  3.08  4.10 

AGE  X  10~5 


C 


1 


PROBABILITY 


APPENDIX  C  GRAPHICS  GENERATED  BY  THE  MODEL 
FIGURE  5.  HISTOGRAM  OF  FAILURE  TIME  FOR  BRIDGE 


r  ■ - a -  “I  ■  1  i  mi, 
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GRAPHICS  GENERATED  BY  THE  MODEL 


FIGURE  6.  HAZARD  RATES 


HAZARD  RATE  = 


#  FAILURES  Ct,t  +  At,) 
**  SURVIORS  AT  t  +  At 
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FIGURE  7.  TOTAL  TIME  ON  TEST  PLOT 


1 


H 

H 
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a 
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c 

CJ 
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V) 


o 
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PROPORTION  OF  FAILURES 


TTTj.  =  TTT^.  _± 


(N  -  I  +  1)  A  Tj 
SCALE 


N 

WHERE,  SCALE  =  ^  (N  -  J  +  1)  A  Tj  AND  I  =  1,  2,  3,  .  .  .  ,  N 

J  =  i 
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appendix  c  graphics  generated  by  the  model 


FIGURE  8.  EXPECTED  COST  OF  REPLACEMENT  POLICY  PER  UNIT  OF  AGE 

a  4,05i 


3.04 


2.03 


1 . 01  il 


Q  ,  QQ  j - 1 - 1 - ! - 

0.02  1.04  2.06  3.08  4.10 

AGE  X  1G~5 


COST  PER  UNIT  OF  AGE  = 


EC COST). 


WHERE,  ECCOST)x  is  the  expected  cost  per  cycle  for  replacement 
at  AGE^.  and  ECCYCLE)^.  is  the  expected  cycle  length  for 

replacement  at  AGE  . 
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DATA  RECORD  FOR  BLOCKS 


TABLE  5.  DATA  ELEMENTS  FOR  BLOCK  RECORDS 


ELEMENT 

TYPE  OF  DATA  ELEMENT 

NAME 

STRING  OF  10  CHARACTERS 

AGE 

REAL  NUMBER 

FAILURETIME 

REAL  NUMBER 

SURVIES 

REAL  NUMBER 

NEXT 

POINTER  TO  A  BLOCK  RECORD 

ARC 

ARRAY  OF  10  POINTERS  TO  BLOCK  RECORDS 

TYPEOFDIST 

USER  DEFINED  DATA  TYPE  C EXPONENTIAL, 

FIXED. 

UNIFORM,  VEIBULL,  NORMAL, 

LOGNORMAL) 

DISTRIBUTION 

VARIANT  RECORD  DEPENDENT  ON  VALUE  OF 

TYPEOFDIS 

EXPONENTIAL  ►  MEAN  TIME  TO  FAILURE 

:  REAL 

FIXED  ►  CONSTANT 

:  REAL 

UNIFORM  ►  MIN,  MAX 

:  REAL 

VEIBULL  ►  SHAPE,  SCALE,  LOCATION 

:  REAL 

NORMAL  ►  MEAN, 

:  REAL 

STANDARD  DEVIATION 

LOGNORMAL  ►  p,  o 

:  REAL 

Table  5.  contains  the  data  elements,  which  make  up  each  record. 

The  element,  NEXT,  combines  with  a  pointer  variable,  FIRSTRECORD, 
to  form  the  singly  linked  list  which  allows  sequential  access  to 
the  block  records.  The  program  constructs  this  list  as  the 
information  elements  of  the  records  are  input  and  insures  that 
each  record  has  a  unique  element,  NAME.  The  tree,  which 
represents  the  blocks  diagram,  is  then  constructed  by  inputing  the 
element,  NAME,  for  the  records  making  up  the  origin  and 
destination  nodes  of  each  arc.  The  program  searches  the  linked 
list  to  find  the  record  associated  with  the  element,  NAME  for  the 
origin  node  and  then  searches  for  the  record  associated  with  each 
destination  node  and  assigns  a  value  to  the  array,  ARC,  which 
points  to  the  destination  record. 
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1 

1 

1 
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